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Abstract 

Effects of the adsorption-desorption process on the immittance response of an electrolytic cell 
are theoretically investigated in the framework of the diffusional Poisson-Nernst-Planck (PNP) 
continuum model, when the generation and recombination of ions is taken into account. The 
analysis is carried out by searching solutions for the drift -diffusion equation coupled to the Poisson's 
equation relating the effective electric field to the net charge density. The effect of different ion 
mobilities on the immittance, i.e., situations with equal and different diffusion coefficients for 
positive and negative ions, are considered. A general exact expression for the admittance in the 
context of the linear approximation is obtained. 
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I. INTRODUCTION 



An important mechanism to explain the role of the ions on the electrical impedance re- 
sponse or immittance response of an insulated material (solid, liquid or gel) may be the 
selective adsorption of ions at the solid electrodes limiting the medium |2|. Thus, to 
build a more complete model to evaluate the electrical impedance of this system it may 
also be necessary to take into account the dissociation of neutralparticles in ionic products, 
the recombination of these ions giving rise to neutral particles [3|-[l5|, and the adsorption- 



desorption process occurring at the interfaces 
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-[l8|. In a recent paper 19|, the importance 



of the dissociation-association process in the determination of the electrical impedance was 
analyzed in detail, but without taking into account the adsorption phenomena, in the frame- 
work of the Poisson-Nernst-Planck (PNP) model. A conspicuous effect of this dissociation- 
association process is the appearance of a new plateau in the real part of the electrical 
impedance of the cell in the dc limit. 

In this paper, we investigate the importance of the dissociation-association phenomenon 
for the electrical impedance of an insulating medium containing ions [20] when the 
adsorption-desorption process is taken into account. Specific adsorption effects at the elec- 
trodes have been already considered by adopting Chang- Jaffe boundary conditions when 



the dissociation-association phenomenon was present [13|,|1J]. The adsorption phenomenon 
considered here is governed by a kinetic equation describing a chemical reaction of first kind 
at the interface liquid medium - electrodes, by imposing the conservation of the mumber 



of particles in entire bulk plus electrodes system [21|, |22|. The mathematical problem to be 
solved is then represented by a set of fundamental equations that are different from the ones 
treated in Ref. jl4|. The strategy adopted here firstly analyze the situation for which the 
positive and negative ions have identical diffusion coefficients, i.e., equal mobilities. This 
scenario may be useful for describing the diffusion of mobile charges in isotropic liquids in 
general. Subsequently, the fundamental equations are faced by means of analytical methods, 
and the solutions are found for the case in which the diffusion coefficients of the positive 
and neutral particles are zero, a scenario eventually more appropriated to the description 
of e.g., an insulating gel. The results we have obtained here are however quite general. In 
fact, they have been obtained in the framework of the PNP continuum model, in which the 
diffusion equation is solved together with the Poisson's equation in order to establish the 
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correct spatial profile of the electrical potential inside the sample. 

The outline of the paper is as follows. In Sec. II, we present the fundamental equations 
of the problem relevant to the drift-diffusion problem for the ions in an insulating medium 
(solid, liquid, or gel), in the presence of the generation-recombination phenomenon, described 
as a first order chemical reaction, and incorporating a kinetic equation at the interface. 
The case in which the diffusion coefficients of the positive and negative ions are equal is 
considered in Sec. III. The particular case in which the positive and neutral particles are 
stuck on the polymer chains forming the gel is discussed in Sec. Ill, where it is shown that the 
adsorption-desorption phenomenon is responsible for a low frequency plateau. In addition, 
the generation-recombination process may induce a plateau for intermediated frequencies 
depending on the thickness kt and on the Debye's length A considered. Some concluding 
remarks are drawn in Sec. IV. 



II. FUNDAMENTAL EQUATIONS 



The fundamental equations established in this Section are relevant to a sample of an 
insulating medium limited by two identical plane-parallel electrodes, a distance d apart. 
The medium, of dielectric constant e, contains impurities that can decompose according to 
the chemical reaction A + C~ , where A indicates the neutral specie, and and 

C~ the positive and negative (monovalent) ions created in the decomposition of A. The 
dissociation constant is indicated by k^, and the association constant by ka- The geometry 
of the sample is such that the electrodes are located at the positions z = ±d/2 of an axis that 
is normal to the surfaces of these electrodes. In addition to the situation treated in Ref . 19| , 
where the sample has the same geometry, we incorporate the adsorption-desorption process 
on the surface of the electrode. 

If we indicate by Nn, Np and Nm the bulk densities of neutral, positive, and negative 
particles, respectively, in the presence of an electric field, of electrical potential V, the 
respective bulk densities of current of the particles are 



oz 



Jp 



d 



dz kuT "dz 



' nSv] 



and 



(1) 
(2) 
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OZ 



keT' '"dz 



(3) 



in which the corresponding diffusion coefficients are Vp and Vm- The equations of 
continuity, stating the conservation of the particles, are 



d 



d 
d 

-Q^jp + ^dNn - kaNpNm, and 
d 

^d.Nn kaNpNm- 



(4) 
(5) 
(6) 



The remaining equation of the model is the equation of Poisson, relating the effective electric 
field to the net charge density. It is obtained by considering that the electrical displacement 
is given by V ■ D = p/e, where p is the bulk density of the charges in the medium. This 
quantity is given by p{z,t) = q{Np — Nm), where q is the ion electrical charge. In the 
one- dimensional case we are considering, the equation of Poisson can be simply written as 



dz- 



■V{z,t) 



{Np - Nm) 



(7) 



The set of fundamental equations 



have to be solved with the boundary conditions 



3n{±d/2,t) = 0, 



d 



Jaiz,t) 

V{±d/2,t) = Vo{±d/2,t), 



and 



(8) 
(9) 
(10) 

r""'^' ^^^^ 
related to the assumption of the adsorption-desorption process and the presence of the ex- 
ternal voltage Vo(t), where a = p refers to the positive ion and a = m to the negative ion. 
As already mentioned, a similar problem was worked out in a different way, by considering 
the Ch ang - Jaffe boundary conditions to take into account specific adsorption at the elec- 
trodes [isl . [l^ . However, it is worth mentionin g ag ain that the equations to be solved here 



are different from the ones considered in Ref. 



14| due to the kind of particles considered 
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in both approaches. In the kinetic equation Eq. ( !TT|) . k and r are parameters describing 
the adsorption phenomenon. This equation simply states that the time variation of the 
surface density of adsorbed particles aa{t) (i.e., of both signs) depends on the bulk density 
of particles just in front of the adsorbing surface, and on the surface density of particles 
already adsorbed. The parameter r has obviously the dimension of time, whereas k has the 
dimension of a length/time. Consequently, if the adsorption phenomenon is present, from 
the kinetic equation, it follows that there is an intrinsic thickness kt. As a final remark, we 
notice that by means of the Eqs. ([4]), (|5]), and ([6]) the generation and recombination of the 
ions are taken into account, whereas by means of the Eqs. ([9]) and ( |TT|) . the adsorption- 
desorption process is incorporated to the dynamics of the distribution of the ions in the 
cell. 

Following the developments performed in Refs. 
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19[, we write = A/'„ + n„, Np 



77 + Up, and Nm = 77 + rim- We assume that the variations in the bulk density of ions due to 
the action of the external field are very small with respect to the values in thermodynamical 
equilibrium. This linear approximation implies that the external difference of potential 
applied to the cell is such that Afn ^ ^n, 77 Up and 77 ^ rim- Likewise, for the 
adsorbed quantities (Wp and ct^), we have Wp = a + Sp and = a + s^, with a ^ Sp 
and a > Sm- Note that Nq = N + Nn and kdNn = kj7 , with 77 = J\f/{1 + 2KT/d) and 
a = [{KT/d)/{l + 2KT/d)]j7d. By accomplishing these approximations, the condition stating 
the conservation of the number of particles takes the form 



[sp + .^) + l'^' ( + ^^^^^ ]dz = 0. (12) 

-d/2 



We also consider the applied potential given by Vo{±d/2,t) = ±(Vo/2) exp{iut), where 
Vq is the amplitude and u the circular frequency of the applied voltage. In the limit of 
small applied voltage (linear limit) the problem was solved in the absence of adsorption- 



desorption process 



found 



19| and the solution without generation and recombination of ions was 



16l-ll8|. In the next section, we obtain analytical solutions for the problem in some 



particular situations, in order to get closed expressions for the impedance and admittance 
of the system. 
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III. SOLUTIONS, IMPEDANCE AND ADMITTANCE 



In general terms, it is possible to look for solutions of the problem in the form 

na{z,t) = rja{z)e^'^^ , with a = n,p,m, 

Sa(t) = aaC^'^^ and 
V{z,t) = (f){z)e''^'. (13) 

Using the approximation discussed above and the ansatz (|T3l) . the fundamental equations 
of the problem can be rewritten as 



iur]n = Vn-j^Vn- kdT]n + kaAf{r]p + r]m) (14) 

- kaAf{T]p + r]ra) (15) 

dz'^ kbi dz"^ 

- kaM{7]p + rira) (16) 

—0 = 

dz"^ e 



-iVp-Vm)- (17) 



Equations ( 1T4|) -( !T7|1 have to be solved with the boundary conditions, at z = ±d/2, 



'^n^Vn = (18) 

-""^di^- - k^^'^Tz'^ = 

^ d qAf ^ d , . 

-T^m-rVm + -r^T^m-r9 = ^^^CTm (20) 

dz kbi dz 

0(±rf/2) = ±1^0/2. (21) 

with cr„ = nrlil + iuJT)!]^. The solution of the problem in the absence of adsorption- 
desorption process has been presented in |20|. In the following, we discuss the particular 
cases where = Vm = T> with P„ ^ V, and Vp = P„ = with Vm = T> ^ 0. The first case 
is rather simple, and could be used to describe systems like a water solution of KCl, close 
to the saturation. The second case could be of some importance in describing the behavior 



Dute to the conduction mechanism, 



23|. 



of gels doped with salt because only negative ions contri 
since the positive ones are stuck on the polymer chains 

Let us analyze the first situation which corresponds to all diffusion coefficients different 
from zero. By substituting in Eqs. (HM - ffT7|) Vp = = T) and P,„ ^ V, it is possible, 
after some calculation, to obtain the equation 



r 



subjected to the boundary condition, at 2; = ±d/2 



V 



ij' +\2q 



(22) 



(23) 



kBT^ 

where ip = rjp — rjm and cr^ = dp — am, with a'^ = ht/{1 + iijjT)ip. These equations lead 
us to obtain the impedance presented in Ref. [21]. This feature implies the generation and 
recombination of ions has no effect on the impedance when the mobility of the positive ions 
is equal to the mobility of the negative ions. Similar result was found in Ref. 20(| in absence 
of adsorption-desorption process. 



Now, we address our attention to the case Dr 



0, with Vm = V T^Q. For this 



case, the set of fundamental equations of the model, namely Eqs. dl]) - ©, becomes 



iurin = -kdVn + kaJfiVp + Vm), (24) 

iUT]p = kdT]n - kj\f{7]p + 7]m), (25) 

- kjj{rip + rim), and (26) 

^0 = -'^{Vp-Vm)- (27) 

The boundary conditions on rjn and rjp are identically satisfied because Vn = Vp = 

imply that the bulk densities of currents for the two types of particles vanish. The remaining 
boundary conditions, taking into account the surface adsorption, are 



T>—rim - --^V—(t) = ±iujam and 
az ksi dz 

<P{±d/2) = ±1^0/2, (28) 
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with 



(^m = JT—. r??m, (29) 

(1 + lUT) 

at z = ±(i/2. From Eqs. and (l25l) . it follows that rjn + rjp = 0, and 

(fcd + A/ fca + ^W) 

for the spatial parts of the variations of the bulk densities of the neutral and positive particles. 
Equations ( l26l) and ( 1271) now become 



and 



where 



^-^ - fG.™ (32) 



zcu + + kaJ\ 



respectively. The solutions, taking into account the boundary conditions (!28|) . are 



Asinh(/3;z) 



(34) 



and 



Gr]m{z) + C z, 



(35) 



with 



. w G iujkj\f 



P = «— + =2 + ; 7^rr\ — • \^^) 



Boundary conditions ( 128|) yield 
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with 



_ VoqPM/ i2ksT) 

u sinh (/3rf/2) + [d£/ (2/3)] cosh {pd/2) ^ ' 

^ Vo {£/(5) cosh {[5d/2) + V sinh {(5d/2) 
2 u sinh (/3rf/2) + [rf^/ (2/3)] cosh (/3d/2) ' 



/3A 2 (1 + zwr) I? 



and A = 2A^, where A = \J eksT / {'iN'q^) is the Debye's screening length. 
From the previous results it is possible to obtain the electric field as follows 

E{z,t) = -§^Viz,t) = -j-^me'^' (40) 

and, consequently, from the Coulomb theorem: E{d/2,t) = — (S(t) + q'cr(t)) /e, where S is 
the surface density of charge on the electrode at z = d/2 and qcr = q{ap — am) is the net 
adsorbed charge at z = d/2. These quantities are relevant to obtain the admittance Y of 
the system. The current is determined by the equation / = SdH/dt, which lead us to the 
result 



, W-lGoosh^/Jjj+^smh^/Jjj+^Cj. (41) 

from which the admittance, Y = I/V, of the sample (cell) may be determined. It is given 

by 



y = Ig cosh ( + sinh (A + ^C]. (42) 

/3Vo \ y 2j l + icoT y 2j qA j ^ ^ 

From this equation it is possible to find the impedance related to this system with 
adsorption-desorption in presence of generation and recombination of ions taking the re- 
lation Z = 1/Y into account. 

Figure [U shows the behavior of the impedance for some significant situations. The black 
(solid) line is the case without adsorption-desorption in absence of dissociation and asso- 
ciation, i.e., K = 0, fed = 0, and ka = 0. The red line is the case k^ = 0.1 with 
ka = 2 X 10^^*^ m^s^^ in absence of the adsorption-desorption. The green (dotted) line 
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iog(/) 

FIG. 1: Real part of the electrical impedance of the cell, R = Re [Z], versus the frequency of the 
applied voltage, / = u)/{2-k), for different values of k, r, k^, and ka- 

corresponds to the case kd = 0.1 s~^, with ka = 2 x 10~^° m^s~^ in presence of the 
adsorption-desorption process, with k = lO^^^ms^^ and r = 50s. The blue (dashed-dotted) 
line corresponds to the case k^ = 0.1 s~^, with ka = 2 x 10~^° m^s~^, k, = 10~^ms~^ and 
r = 50s. One observes that the green dotted and the blue dashed-dotted lines take the effect 
of adsorption-desorption and association and dissociation of ions into account. They evi- 
dence that the adsorption-desorption process may play an important role at low frequency 
while the generation and recombination of ions has predominant effect for intermediated 
range of frequencies. 

Figure [2] shows the behavior of the imaginary part of the the impedance for the same set 
of parameters used in Fig. 1. Note that the imaginary part of the impedance is not sensible 
to the adsorption-desorption when < A, as shown by the green dotted line. On the other 
hand, for kt > X the effect of this process is evidenced (see the blue dashed-dotted line). 
These features can be established by analyzing the low frequency limit of the admittance. 
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iog(/) 



FIG. 2: Imaginary part of the electrical impedance of the cell, X = — Im [Z], versus frequency of 
the applied voltage, / = uj/{2-ir), for different values of k, t, k^, and ka (see the text). 

Indeed, in the a; — )■ limit, the asymptotic behavior of the admittance is given by 



Y ^ iuY- {iuYV, (43) 

where 



F ~ S 
A 



KT H I d H 



and 



(44) 



Y 



Se 



AV (XF)' 



[dkJ^V - dX^FH 



+ 2ktF^ 



T-2 



{2V + dK)r- AX 



+ 2X 



X 



kaVj7 + FH 2dKT - 3A 



coth 



d H 



2X\ F 
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FIG. 3: Parametric plot of X versus R for different values of k, t, k^, and ka- 

+ d (SFHX' - kaMv) coth (^i^f j I . (45) 

with F = kd + kj\f and H = kii + 2kJ\f. By using Eq. (143|) . it is possible to show that the 
asymptotic behavior of the impedance for low frequency is given by i? = Re (Z) Y /y'^ 
and X = Im [Z] ^1/ {iuY^ . In Fig. [3l a parametric plot of X versus R is reported. As the 
frequency increases, the effect of the adsorption-desorption is evidenced before the process 
of generation and recombination of ions and it may also be present in the high frequency 
limit. In Figs. |l]and[5l we illustrate the behavior of e' and e" for different situations with 
adsorption-desorption and generation and recombination of ions. In particular, the set of 
parameters used to plot these figures is the same of Fig. [TJ In all these figures, we are 
also considering the ions monovalent with q = 1.6 x 10~^^ As and the temperature is such 
KbT/q = 0.025^. The geometric parameters of the cell are supposed to be = 25yum and 
S = 2 X 10~^m^. We also assume that the dielectric constant of the insulating material is 
e = Q.Teq, the bulk density of impurities in thermodynamic equilibrium is Ao = lO^^m"^, 
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FIG. 4: e'/eo versus f = uj/{2-k) for different values of k, r, kci, and ka- 



and the diffusion coefficient of the negative ions is P = 8.2 x 10 ^^m^s ^ 



IV. CONCLUDING REMARKS 



We have investigated the effect of the adsorption-desorption process on the electrical 
response of an electrolytic cell (a finite-length situation) when the generation and recombi- 
nation of ions are present. We have performed this analysis by considering two situations. 
The first of them is characterized by equal mobilities for positive and negative ions, which 
shows that the electrical response is not sensible to the generation and recombination of ions. 
It is only influenced by the adsorption-desorption process. The other case considers different 
mobilities for the ions. It is shown that the electrical response can be influenced by both 
effects, i.e., adsorption-desorption and generation and recombination of ions. In this sense, 
one concludes from Fig. [1] and Fig. [2] that the system may exhibit different behaviors for 
the electrical response. In fact, in Fig. [T]the generation and recombination process governs 
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FIG. 5: e"/eo versus / = ti;/(27r) for different values of k, r, kd, and ka- 



the behavior of the system for intermediate frequencies. Indeed, one can see the presence 
of initial part of a second plateau and, after it, the third plateau due to the adsorption- 
desorption process of ions (the green dotted line). The blue dashed-dotted line is essentially 
dominated by the adsorption-desorption process. In Fig. [21 it is shown that the imaginary 
part of the impedance may not manifest the effect of the adsorption-desorption process as 
in the real part depending on the thickness kt and the Debye length A (see, for example, 
the green dotted and blue dashed-dotted lines). On the other hand, this effect and the 
generation-recombination of ions are evidenced for the blue dashed-dotted line which, for 
the real part, is essentially governed by the adsorption-desorption process. Thus, the values 
of the thickness kt and the Debye length A are useful to inform us about the influence of the 
surface and bulk effects on the dynamic of the ions in an electrolytic cell. We hope that the 
present theoretical approach, due to its generality, may be useful to discuss the immittance 
response of an electrolytic cell when both, the process of generation and recombination of 
ions and the adsorption-desorption phenomenon, are taken into account by means of a set 
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of fundamental equations charactering a continuum diffusional model. 
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